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ABSTRACT 

Non-thermal hard X-ray and high-energy (HE; 1 MeV < E < 100 GeV) 
7-ray emission in the direction of r] Carina has been recently detected using 
the INTEGRAL, AGILE and Fermi satellites. So far this emission has been 
interpreted in the framework of particle acceleration in the colliding wind region 
between the two massive stars. However, the existence of a very fast moving 
blast wave which originates in the historical 1843 "Great Eruption" provides an 
alternative particle acceleration site in this system. Here we explore an alternate 
scenario and find that inverse Compton emission from electrons accelerated in the 
blast wave can naturally explain both the flux and spectral shape of the measured 
hard X-ray and HE 7-ray emission. This scenario is further supported by the 
lack of significant variability in the INTEGRAL and Fermi measured fluxes. 

Subject headings: stars: individual(Eta Carina) — acceleration of particles — 
shock waves — radiation mechanisms: non-thermal 
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1. Introduction 

For a long time t] Carina - one of the most peculiar objects in the Milky Way - was 
believed to be a hypergiant luminous blue variable (LBV) star. Recent observations, however, 
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Tj Carina experienced a historical outburst in the 19**^ 
of gas which moves outwards at an average speed of ~ 650 km s~^ 
This material formed the Homunculus Nebula which is expanding with 
i (4 — 10) X 10^^ erg - roughly 10% of the energy released by a supernova 



The optical flDaminelilll996l : lOamineh et al.ll2000f ) and IR flWhitelock et al.lll994 \2004 \ 
lightcurves suggest a binary period of 5.5 years, r noreover, the analysis of the X-ray lightcurve s 
points to an highly eccentric orbit of e ~ 0.9 (ICorcoran et al.l l200ll : lOkazaki et al.ll2008l ). 
rj Carina is a colliding wind binary (CWB) with powerful winds produced by both mem- 
bers of the systern . The primary has a very high mass-loss rate of Mi > 5 x 10~^ Mq yr~^ 
(IHillier et al.ll200ll : iParkin et al.ll2009l ) and a terminal wind velocity of vi ~ (500—700) km s^^ 
whereas the secondar y has a fast low-density wind with V2 ~ 3000 km s^^ and M2 ~ 
1.0 X 10~^ Mq yr~^ (jPittard fc CorcoranI |2002| ) . When supersonic expanding winds col- 
lide, they form a wind-wind interaction zone where charged particles can in principle be 
accelerated to high energies via the di ffusive shock acceleration (DSA) process (see e.g. 
Eichler &: Usovlll993l : lReimer et al.ll2006l ). Observationally, particle acceleration in the wind- 
collision region of massive binary systems such as 77 Carina is suggest ed by the detection of 
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non-thermal radio emission in several objects (see 
indicates the existence of high-energy (HE; 1 MeV < E < 100 GeV) electrons in these 
systems. 

There is also evidence for the presence of relativistic particles in ry Carina with non- 
thermal X-ray emission recently reported by the INTEGRAL collaboration. The authors 
obtain an average l uminosity of 7 x 10^^ ergs s~^ in the hard X-ray band (22 — 100) keV 



(iLeyder et al.ll2008l ). In the HE domain, a source spatially coincident with the y Carina 



posit ion, and possibly variable, was reported by the AGILE collaboration (ITavani et al. 
20091 ) . A steady source in the same direction (including systematical and statistical errors 



see Fig. [T]and Tab. [T]) was confirmed using the LAT instrument aboard the Fermi satellite 



^ Gomez et al.l |20ld ) conclude from submillimetre observations that more than 40 Mq could have been 
ejected in the giant outburst. We restrict ourselves to the more conservative estimate of the kinetic energy 
in the expanding material. 
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flAbdo et alJl2009l . l2010ah . The luminosit y in the HE domain between 100 MeV and 100 GeV 
as obtained from the energy flux given in lAbdo et al.l (l2010al ) for a source distance of 2.3 kpc 
is 1.6 X 10^^ ergs s~^, a factor of 20 higher than the INTEGRAL luminosity. 

Given the large population of HE 7-ray sources detected by Fermi, with a significant 
fraction associated with pulsars, the chance probability of source confusion with a pulsar 
or another HE 7-ray source must be examined. Within a circle of 5° radius, 19 ATNF 
pulsars with a spin-down luminosity E > 3 x iq33 q^^/s ( Manchester et al.l 120051 ) - the 



sensitivity limit of t he Fermi 7-ray pulsar catalog (lAbdo et al. 



2010bl ) - and 12 Fermi sources 



( lAbdo et al.ll2010al ) have been detected. Thus, the probability to find any of the HE 7-ray 
sources or ATNF pulsars within the Fermi error circle at the rj Carina position is estimated 
to less than 0.2%. Moreover, whereas pulsars typically have a cut-off in the energy spectrum 
< 10 GeV, the spectrum of the Fermi source positionally coincident with rj Carina shows 
significant flux above 20 GeV, making an association with a pulsar unlikely. 

The lightcurves from all mentioned instruments, together with lower-energy X-ray data 
from RXTE, which probes the colliding wind region (CWR), are shown in Fig. [2J Note, 
that the flare reported by AGILE around phase p ~ 0.96 was not confirmed using the more 
sensitive Fermi-LAT . Although INTEGRAL detected no significant hard X-ray emission 
during periastron passage, the upper limit is consistent with steady emission. The absence 
of significant variability in the 50 keV and GeV regime as shown in Fig. [2] is surprising in a 
CWB picture, especially during periastron passage where a collapse of the CWR is expected 
and hence no particle acceleration should occur. 

The detections of hard X-ray and HE 7-ray emission in rj Carina have so far been 
interpreted in the framework of particle acceleration in the colliding winds of the two massive 
stars. However, the CWR of 77 Carina is not the only place within the Fermi, AGILE 
and INTEGRAL error circles, where particles can be accelerated via DSA to high energies. 
Recent observations show that rj Carina is surrounded at a distance of ~0.25 pc by a very 



Instrument 


Energ 


y range 


Position (J2000) 
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INTEGRAL 


22 keV 


- 100 keV 


10^45°'02^ -59°43'38" 


2.8 


Fermi 


30 MeV 


- 30 GeV 


10M5™13^ -59°42'21" 


3.6 


AGILE 


100 MeV 


- 100 GeV 


10h44™49^ -59°44'54" 


24 



Table 1: Best fit positions of the sources detected in the hard X-ray and HE 7-ray band. 
The 90% (stat.) error on the position of the INTEGRAL source as well as the 95% (stat. -|- 
sys.) uncertainty on the source position of Fermi and AGILE are given. 
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Fig. 1. — Color-coded image of rj Carina with X-ray emission (0.5 — 11 keV) as measured 
by Chandra (ObsID: 1249) in linear scale shown in blue and optical emission as measured 
with the HST in log scale in red. Overlaid in red and green are the 2.8 ^ INTE GRAL (90% 



stat.- lLevder et all (120081 )) and 3.6' LAT (95% stat. + svs.. lAbdo et all fl2010af )) positional 
uncertainties of the hard X-ray and HE 7-ray source positions, respectively. The positional 
uncertainty the AGILE source encompasses the entire field. The best fit positions and 
corresponding errors are summarised in Table [H 
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Fig. 2. — rj Carina lightcurve as a function of orbital phase as measured by RXTE (black), 
INTEGRAL (blue), AGILE (green) and the Fermi-LAT instrument (red). 
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fast moving blast wave. The giant outburst of 1843 (also known as the "Great Eruption") 
produced ma terial which is now moving ahead of the Homunculus Nebula at speeds of 3500 — 
6000 km s~^ (ISmithll2008l ). This blast wave currently overruns the "Outer Ejecta" a ring-like 
structure of material which orig inates from an ejection of mass from rj Carina ~ 500 — 1000 
years ago (IWalborn et al.l 119781 ) . 



The existence of the fast-moving material doubles the estimate of the total kinetic 
energ y of the giant outburst and mimics a low-energy supernova remnant (SNR) shell ([Smith 
20081 ) with a blast w ave moving into t he ISM with v elocities com parable to the historical 
supernovae RCW 86 Jvink et al.l[2006h and SN 1006 Jvinkll2005h . Fig. □ shows the X-ray 
emission as measured with Chandra and optical emission measured with the Hubble Space 
Telescope (HST). The 95% (stat. + sys.) and 90% (stat.) error circles of Fermi-LAT and 
INTEGRAL, respectively, are overlaid. They enclose the rj Carina position. In this work the 
possibility of particle acceleration in the expanding blast wave of 1843 's Great Eruption is 
examined and leptonic and hadronic origins of the observed 7-ray emission are investigated 
in detail. 



2. Particle acceleration in the blast wave 



SNRs start their evolution by a free-expansion phase of the ejecta and enter the pressure- 
driven Sedov phase, once the mass of displaced material roughly equals the ejecta mass. The 
transition between the two phases in the remnant surrounding t] Carina is expected to occur 
after 290 (500) njj^^^ years depending on the ambient density tih, for a blast wave speed of 
6000 (3500) km/s, respectively. The density of the medium surrounding the blast wave can 
be estimated by adopting two steady stellar winds which filled the volume around r] Carina 
prior to the Great Eruption. For typical wind speeds of 500 km/ s and a typical mass-loss 
rate of 5 - 10 x 10"^ M© yr'^ JmUier et al.l I2OO1I : IParkin et all mOS) the density of the 
circumstellar medium would be of the order of 50 cm~^ at a distance of 0.25 pc. Hence, the 
blast wave is currently in the transition between free-expansion and Sedov phase, implying 
that efficient particle acceleration is possible. 

The dist ance of 0.25 p c of the blast wave from 77 Carina used here is motivated by the 
conclusion of ISmithI (120081 ) that parts of the Ou ter Ejecta are curren tly being overrun by 
the blast wave and by X-ray observations, where ISeward et al.l (1200 ll ) used a ring with the 
same radius to describe the soft X-ray shell coincident with the Outer Ejecta. On the other 
hand, for a steady shock speed of 3500 km/s and an age of 167 years, the expected distance 
of the blast wave to rj Carina is 0.6 pc. A wind solution of the fast moving material, which 
would result in an increasing velocity with time, could explain such a difference. To account 
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for the different distance estimates, the spectral energy distributions (SEDs) for both values 
are shown in Section |H 

Relativistic electrons and protons will gain energy when crossing the shock front of 
rj Carina's blast wave. In the limit of a strong shock, the particle crossing time At is 
determined by the diffusion coefficient k, expressed in terms of the Bohm diffusion coefficient: 

1 n23 2 -1 / -E'ceV 

K = riKBohm = 10 T] Cm S ' 



the shock speed Vs, the particle energy EceV and the magnetic field in the shock region B^q 
as: 

At = ^^^^^ = 1.75 X 10^ s ( ^ ^ . (2) 

CVs V^',,l03 km s-^B^G J 

The energy gain per shock crossing is given by AE/E = Vg/c and, hence, the acceleration 
timescale race = E / {AE / At) ^ rjn/vl: 

« 5 X 10' = r, (^) . (3) 

No magnetic field measurements are available in the region of the blast wave. However, for 
the Homunculus Nebula different estimates exist and are used in the follow i ng as a rough 



guide for the magnetic field strength in the region of interest. lAitken et al.l ( 119951 ) derived 
from polarization measurements, which were based on dust grain alignment, magnetic field 
strengths of C?(/iG-mG) in the Homunculus Nebula, depending on the underlying process. 
However, the dust shell of the Homunculus appears to be neutral and is likely composed 



of silicates and Fe rather than typical ISM dust grains (IGail et al.l l2005l ). Hence, for the 
purpose of this work, a lower magnetic field strength of 10 /xG is assumed. In the Bohm limit 
and without energy losses, it would take about a year to accelerate an electron or proton to 
an energy of 100 GeV in a 10 /xG field given a shock speed of 3500 km s~^. This acceleration 
timescale is short compared to the time since the giant outburst (167 years) and implies that 
charged particles can potentially be accelerated to TeV energies. 



3. Origin of the 7-ray emission 



Protons and heavier nuclei accelerated in the expanding blast wave will interact with 
the ambient medium, producing 7r*^-decay 7 rays above ~300 MeV. The timescale for inelas- 
tic proton-proton interactions is somewhat energy- dependent but can be approximated as 
(jAharonianll2004 ): 



Tpp = 3 X 10^ 
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yrs 



(4) 
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where uh is the density of hydrogen at oms per cm~ ^ . From a veil of [O III] emitting material 
surrounding rj Carina and its ejec ta ISmith et al.l ( 120051 ) estimate a mean density of 400 
particles per cm~^ and lSmithI (120081 ) finds 500 cm~^ in this region, implying = 6 — 8 x 10^ 
years. This timescale is much longer than the age of the remnant and hence the p-p channel 
is radiatively inefficient. To explain the observed level of 7-ray emission with this density of 
target material an energy of about L^e^Tpp ~ 3 x 10"^^ erg in hadronic cosmic rays (e^ ~ 0.1 
is the fraction of the proton energy released in 7 rays in a typical interaction) would be 
required. This is roughly 5% of the kinetic energy of the blast wave, quite plausible within 
the framework of DSA. However, it is hard to explain the Fermi emission below 300 MeV, and 
especially the INTEGRAL emission, in this scenario. Furthermore the required maximum 
energy of protons as indicated by the curvature in the Fermi spectrum is well below the 
expected maximum acceleration energy associated with either the age or size of the system. 

Contrary to hadronic cosmic rays, accelerated electrons will encounter significant en- 
ergy losses during the acceleration process. They will predominantly lose energy via in- 
verse Compton (IC) scattering on radiation fields present in the acceleration region, via 
synchrotron radiation in magnetic fields, and by Bremsstrahlung in interactions with am- 
bient material. The IC cooling time in the presence of a black-body target radiation field 
with energy density Urad . ev cm-3 for an electron of energy Eq^v is (IModerski et al.l 12005 
Hinton fc HofmannI 120091 ) : 



rioss ~ 3.1 X 10« yrs Ur' 



rad, eV cm~3 GeV 



/kn 



KN • 



(5) 



Assuming that the radiation field of the whole Homunculus Nebula is concentrated in a point 
at 0.25 PC distance to t he blast wave and follows a black-body spectrum of 260 K temperature 
(jCehrz fc Smithlll999l ). the energy density in IR photons which serve as target radiation field 
is f/rad = 8 X 10^ cV cm~^, several orders of magnitude higher than the density of the cosmic 
microwave background (0.27 eV cm~'^) or typical ISM photon fields. For infrared target 
photons the Klein-Nishina (KN) suppression of the IC cross-section for electron energies 
E < 100 GeV is a less than 20% effect and is not considered in the following. The relevant 
IC cooling time for electrons in the blast wave region is therefore rioss ~ 4000 yrs E^ly. In 
a 10 /iG field, the synchrotron-loss timescale for 100 GeV electrons is of O(IO^) years: IC 
losses will dominate unless B > 2mG. In dense media, Bremsstrahlung may also contribute 
significantly to the energy-loss rate of electrons. The estimated mean density of 400 — 500 
particles per cm~^ implies a typical Bremsstrahlung loss time of TBrcms ~ 8 — 11 x 10^ years. 
Therefore, at > GeV energies, IC emission is very likely to be the dominant emission process 
for electrons in this system. 



The maximum electron energy is determined (for cooling-limited acceleration) by the 
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balance between radiative energy losses and the rate of energy gain via DSA, when rioss = Ta 



^e,ma. ~ 50 GeV V^,^^, km s'M / ^ " (6) 



V 



This maximum energy for the measured blast wave speed is consistent with the HE 7- 
ray emission found by Fermi for moderate magnetic fields of 10 and diffusion close to 
the Bohm limit. Moreover, the energy-loss timescale is much shorter than the age of the 
remnant, implying that acceleration of electrons is indeed cooling limited in this system. In 
this scenario most of the energy injected into the acceleration process emerges as IC radiation 
in the GeV domain. For such an efficient process the required energy input in electrons is 
modest: within an order of magnitude of L^tgys ~ lO''^ erg. IC emission is expected to 
dominate over p-p emission in this system as long as E^. > IO^'^-Ecr- 



Vannoni et al.l ( 120091 ) performed numerical calculations where they studied acceleration 
and radiation of electrons in an radiation-dominated environment such as the one investigated 
in this work. Compared to the magnetic field of 10 fiG used here, the authors obtain for 
magnetic fields of 100 fiG maximum particle energies in the TeV regime, where KN effects 
become very significant. 



4. Spectral Energy Distribution 

Figure |3] shows the measured SED of the 1' region around rj Carina together with 
calculated synchtrotron and IC emission for the leptonic scenario outlined above (solid 
lines). The calculation is a single- zone time-dependent numerical model as described in 



Hinton fc AharonianI (120071 ). 



Electrons accelerated via the DSA process are expected to follow a power-law in energy 
dN/dEe ~ E~^^ with index Fg ?5i 2 and produce an IC (and synchrotron) 7-ray spectrum 
with spectral index F^ = — (Fg + l)/2. For continous injection IC cooling in the Thomson 
regime results in a spectral break with AF^ = 1/2. In the case of the Great Eruption 167 
years ago, the IC break energy is expected to lie at ~ 0.4 GeV for a blast wave distance of 
0.25 pc. As Figure |3] shows, for this distance the observed break between the INTEGRAL 
and Fermi spectra occurs somewhat below this energy, close to the expectation for the 
higher radiation field corresponding to a smaller distance of 0.15 pc. This value is broadedly 
consistent with expectations since the relevant radiation field is the average value encountered 
over the history of the object. Furthermore, the system is far from radially symmetric and 
parts of the shell lie at significantly smaller than average distances to the Homunculus. The 
Fermi data indicate a high energy cut-off at an electron energy of ~ 100 GeV, very close to 
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the expected cut-off for IC-cooling-limited acceleration. 

The dashed and dotted-hnes in Figure [3] demonstrate the effect of adjusting the main 
parameter of the model. For the dashed curves the radiation field is decreased to that 
expected at 0.6 pc (rather than 0.25 pc) from t] Carina. In this case an increase in energy 
in electrons by a factor of 2.3 would be required to match the Fermi flux. For the dotted 
curves the radiation field is increased to that expected at 0.15 pc from t] Carina. All curves 
were produced assuming a magnetic field of 10 /iG. Magnetic fields higher than ~ 20 /xG are 
excluded in this scenario by the limit on non-thermal radio emission from Molonglo. For a 
blast wave distance of 0.6 pc, the shape of the resulting model spectrum reproduces the data 
in the HE domain less well. A model with a realistic time and space dependent radiation 
field would be expected to fall in the region defined by the three curves shown. 



5. Summary 

It appears that the properties of the high-energy non-thermal emission of rj Carina are 
consistent with an origin at the high-velocity blast wave of the Great Eruption of 1843. As 
Figure [3] demonstrates, the position of the spectral break between the Fermi and INTEGRAL 
domains, and the cut-off energy measured using Fermi, can be explained as a cooling break 
and an acceleration limit, respectively. 

A prediction of this model is an extended emission region and hence non-variability of 
the signal. The non-variability of the source in the 7-ray regime is supported by the Fermi 
observations. While Fermi and INTEGRAL do not provide sufficient angular resolution to 
resolve the blast wave (see Fig. [T]), high resolution radio observations will be able to resolve 
non-thermal radio emission from this acceleration site. Detection of a radio shell at the 
location of the shock would support the proposed scenario. 

The time elapsed since the ejection of the very fast material by rj Carina is comparable 
to the age of the youngest Galactic SNR G 1.9+0.3 and between the ages of the remnants 
of SN 1987A and CasA. It may provide important insights into particle acceleration and 
production of non-thermal radiation during the early stages of remnant evolution where 
only very few objects can be studied. 
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- 11 - 




lO"*^ 10"" 10"^ 1 10^ 10" 10*^ 10" 10^° 10^^ 10^" 

Energy (eV) 



Fig. 3. — Measured and calculated SED for the region within a few a rcminutes of rj Carina, 
showing the radio flux measured by Molonglo (iMurphy et al.l 120071 ) as an upper limit on 
the synchtrotron flux, therma l X-ray emission as me asured by MEGS, which probes the 
wind-wind iri t eracti on region, (jviotti et al. 2002, 2004), the INTEGRA L measurement from 



Leyder et al.l (120081 ) and the Fermi-L AT detection ( lAbdo et al.ll2010al ). The sensitivity of 



the H.E.S.S. instrument at TeV energies is shown for comparison. The curves show a single 
zone time- dependent model for continous injection of electrons over the 167 year history 
of the system and a maximum acceleration energy of 110 GeV. A radiation density of 8 x 
10^ eV cm~^, appropriate for a region at 0.25 pc from r] Carina, is used for the solid curves, 
decreased to 1.4 x 10"^ eV cm~^, corresponding to 0.6 pc distance, for the dashed curves and 
increased to 2.2 x 10^ eV cm~^, corresponding to 0.15 pc distance, for the dotted curves. A 
magnetic field strength of 10 fiG is assumed. The Bremsstrahlung component for each model 
is depicted in green. 



- 12 - 



German Ministry of Education and Research (BMBF) is acknowledged. 

REFERENCES 

Abdo, A. A., et al., 2009, ApJS, 183, 46 



Abdo, A. A., et al., 2010, http://fermi.gsfc.nasa.gov/ssc/data/access/\lat/lyr_catalog/lFGL_ca' 



Abdo, A. A., et al., 2010, ApJS, 187, 460A 

Aharonian, F. A., 2004, Very high energy cosmic gamma radiation: a crucial window on the 
extreme Universe, World Scientific Publishing, 2004 

Aitken, D. K., et al., 1995, MNRAS, 273, 359A 

Corcoran, M. F., et al., 2001, ApJ, 547, 1034 

Damineh, A., 1996, ApJ, 460, L49 

Damineh, A., et al., 2000, ApJ, 528, LlOl 

De Becker, M., 2007, A&A Rev., 14, 171 

Eichler, D. and Usov, V., 1993, ApJ, 402, 271 

Gail, H.-R, et al., 2005 ASP Conference Series, 332, 317G 

Gehrz, R. D. & Smith, N., 1999 ASP Conference Series, 179, 251 

Gomez, H. L., et al., 2010, MNRAS, 401L, 48G 

Hillier, D. J., et al, 2001, ApJ, 553, 837 

Hinton, J. A. & Aharonian, F. A., 2007, ApJ, 657, 302H 

Hinton, J. A. & Hofmann, W., 2009, ARA&A, 47, 523 

Leyder, J.-C, Walter, R. & Rauw, G., 2008, A&A, 477, L29 

Manchester, R.-N., et al., 2005, AJ, 129, 1993M 

Moderski, R., et al., 2005, MNRAS, 363, 954M 

Murphy, T., et al., 2007, MNRAS, 382, 382 - 392 



Nielson, K. E., et al, 2007, ApJ, 660, 669N Aharonian, F. A., 2007, ApJ, 657, 302H 
Okazaki, A. T., et al., 2008, MNRAS, 388L, 390 
Parkin, E. R., et al., 2009, MNRAS, 394, 1758P 

Pittard, J. M. & Corcoran, M. P., 2002, A&A, 383, 636 2002, A&A, 388, 335 

Reimer, A., Pohl, M. & Reimer, O., 2006, ApJ, 644, 1118 

Seward, F. D., et al., 2001, ApJ, 553, 832 

Smith, N., et al., 2003, AJ, 125, 1458 

Smith, N., Morse, J. A. & Bally, J., AJ, 130, 1778 

Smith, N., 2008, Nature, 455, 201 

Tavani, M., et al, 2009, ApJ, 698, L142 

Vannoni, G., Gabici, S. & Aharonian, F. A. 2009, A&A, 497, 17V Conference Series, 774, 
241 

Vink, J., 2005, American Institute of Physics Conference Series, 774, 241 
Vink, J., et al, 2006, ApJ, 648, L33 
Viotti, R. F., et al., 2002, A&A, 385, 874 
Viotti, R. P., et al., 2004, A&A, 420, 527 

Walborn, N. R., Blanco, B. M. & Thackeray, A. D., 1978, ApJ, 219, 498W 
Whitelock, P. A., et al., 1994, MNRAS, 270, 364 
Whitelock, P. A., et al., 2004, MNRAS, 352, 447 



This preprint was prepared with the A AS macros v5.2. 



